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Abstract 

The question of whether occupation numbers and momentum distributions of nucleons in nuclei 
are observables is considered from an effective field theory perspective. Field redefinitions lead 
to variations that imply the answer is negative, as illustrated in the interacting Fermi gas at low 
density. Implications for the interpretation of (e, e'p) experiments with nuclei are discussed. 
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The advent of new continuous beam electron accelerators (CEBAF, Mainz, MIT-Bates, 
NIKHEF) permits experiments that probe hadronic matter in both inclusive and exclusive 
reactions with unprecedented precision. These experiments are expected to deepen our 
understanding of nuclear structure and reaction mechanisms by measuring the response of 
nuclei to electroweak probes over a wide range of energy and momentum transfers ||T], 
Exclusive electron scattering experiments at large momentum transfer, such as the knockout 
of a proton in (e, e'p) on a nucleus, are particularly important. 

It is often claimed that occupation numbers and/or momentum distributions of nucleons 
in nuclei can be extracted from these experiments. Their extraction from the data is built 
on the extreme impulse approximation but is obscured by the need to consider final state 
interactions and meson exchange currents. The interpretation and comparison with theory 
is made based on a given nuclear Hamiltonian, but from the perspective of the underly- 
ing theory of the strong interaction, QCD, there is no unique or preferred Hamiltonian. 
Rather, there are infinitely many such low-energy effective Hamiltonians that are related by 
field redefinitions that leave observables unchanged (e.g., see Ref. 0]). What happens to 
occupation numbers under such transformations? 

A powerful framework to study low-energy phenomena in a model-independent way is 
given by effective field theory (EFT) 0, H, |^. The underlying idea is to exploit a separation 
of scales in the system. For example, if the typical momenta k are small compared to 
the inverse range of the interaction 1/R, low-energy observables can be described by a 
controlled expansion in kR. All short- distance effects are systematically absorbed into low- 
energy constants using renormalization. The EFT approach allows for accurate calculations 
of low-energy processes and properties with well-defined error estimates. 

In this letter, we explore from an EFT perspective the question of whether occupation 
numbers and momentum distributions of nucleons in nuclei are observable. In an EFT, 
observables are characterized by invariance under local field redefinitions. If a quantity 
depends on the particular representation of the Lagrangian C (beyond the level of truncation 
errors), it is not an observable. Off-shell Green's functions for scattering processes in the 
vacuum, for example, can be changed by field redefinitions. On-shell Green's functions, 
which correspond to S-matrix elements, however, are unchanged |P, 0. 

To address the issue most cleanly, we focus on the question of whether occupation numbers 
in a homogenous medium at finite density are observables in the framework of EFT. In a 
general finite system, occupation numbers and momentum distributions are very different 
quantities. In a homogeneous system, however, they are equivalent. Since there are no 
asymptotic states in an infinite medium, however, the usual analysis for field redefinitions 
does not directly carry over to in-medium observables. While the analysis can be extended to 
thermodynamic observables like the energy density or particle number ||^, this extension is 
not obvious for other quantities, such as the momentum distribution or occupation numbers. 

We use the interacting Fermi gas at low density as a laboratory to illustrate a funda- 
mental and generic problem with the definition of momentum occupation numbers. For a 
given representation of the Hamiltonian, one can define an operator al.a^^ that counts parti- 
cles/holes with momentum k and gives the expected result in the noninteracting limit. As 
discussed below, this operator is not derived from a global symmetry, in contrast to the total 
number operator. In the EFT framework this is a problem, as there is no preferred form of 
the effective Lagrangian. 

To highlight this problem, we consider both the particle number and the momentum 
distribution n{p). We describe the system using a local Lagrangian for a nonrelativistic 
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fermion field with spin independent interactions that is invariant under Gahlean, parity, and 
time-reversal transformations: 
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where V is the Gahlean invariant derivative and H.c. denotes the Hermitian 

conjugate. A convenient and transparent regularization scheme is dimensional regularization 
with minimal subtraction; see, e.g., Ref. for details. In this scheme the coefficients are 
simply Co = Ana/M and C2 = Coare/2 where a is the s-wave scattering length and the 
effective range. 

We generate an infinite, equivalent class of Lagrangians Ca by performing the field redef- 
inition 

^ ^ ^ + ^(^V)^ , - ^^n^V) , (2) 

in C. The factor is introduced to keep the arbitrary parameter a dimensionless. A is 
the breakdown scale of the EFT and the additional factor of Att is introduced for convenience 
§. We obtain for £„: 

= £ - ^2Co(^V)' 
Aira 

+ ■■■ + 0{a') , (3) 

where higher-order two- and three-body terms, all four- and higher-body terms, and terms 
of O(a^) have been omitted. 

For a = we recover the original C. The Lagrangian Ca contains additional vertices, 
including an off-shell vertex, but gives exactly the same energy density and particle number 
as the Lagrangian C. In Ref. ||^ it was illustrated how the necessary cancellations occur 
in general and for this particular example. Furthermore, it was shown how the number 
operator must be constructed as the conserved charge of the Noether current associated 
with the U{1) phase symmetry: 

^{x) — > e-*V(a;) and tp^ix) — > e'^^^x) , (4) 

under which is invariant. One can conveniently identify the Noether current by promoting 
(j) to a function of x and considering infinitesimal transformations with 

Ca ^ Ca[lp,1p^](p{x)] . (5) 

Then the number density operator is given by [§: 

iV"^^^£„[^,^t.^(^)] = ^t^ + ^2(^V)' . (6) 

The particle number itself is given by the spatial integral of A^". In a uniform system, 
however, the difference is simply a factor of the volume and it is convenient to refer to N"' 
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FIG. 1: Diagrams contributing to S*(u;,k) to 0{C^) [(a) and (b)] and at 0{aCo) [(c) to (h)]. 



as the number operator. Note that Eq. (P) differs from the naive expectation ip'^ip for a ^ 0. 
In the appendix of Ref. it was demonstrated how the contributions from the additional 
vertices in Ca and the additional term in A^" cancel order- by-order in a. Consequently, the 
total particle number N is unchanged by field redefinitions, as expected for an observable. 

Matters become more complicated for the momentum distribution n{k) since the corre- 
sponding operator is not simply related to the conserved charge of a Noether current. In 
the literature, the operator 

nk = alay^ , (7) 

where destroys a fermion of momentum k, is universally adopted [|T], 0, [1^, This 
definition gives the correct result in the noninteracting limit, n{k) = g9{kp — /c), where g is 
the spin degeneracy factor. If there were a preferred form of the Lagrangian/Hamiltonian, 
as is usually assumed, Eq. (|^ would uniquely determine the momentum distribution. In 
terms of the field operators 

*W=/(S5»--X. *'W = /(05»-*^M. (8) 

we can write ci\xi]^ = J d^xnk{x.) with 

nk(x) = J dS e'^-^i^^^ + y)V^(x) , (9) 

which is nonlocal in coordinate space. Using the definition of the one-particle Green's 
function fig, 

«(x.,x-.0^ <'^°l^l^'''-,''f''-'''"l"°> , (10) 

with |\E'o) the exact ground state of the system and 

?/Ji^(x,t) = e^^*?/^(x)e-*^*, (11) 
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FIG. 2: Schematic picture of the occupation number as a function of momentum in a uniform 
Fermi system with no interactions (dashed hne) and including leading correction from interactions 
(solid line) (cf. Ref. [|lO|). The square with a cross denotes an insertion of Eq. (jlj). 



a Heisenberg operator, the expectation value of nk(x) can be written as |TT 

n{k) = (nk(x)) = YiM-i)g I ^e-^G(^,k) . (12) 

The momentum distribution n{k) for the hard sphere Fermi gas to second order in kpa has 
been calculated using this definition in Refs. [l^, . 

With the definition in Eq. ([12|) , the explicit expressions for n{k) in Refs. |jl3|, ^ 
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reproduced in the EFT by setting a = and calculating the one-particle Green's function 
to 0{Cq). The one-particle Green's function G'(u;,k) is related to the proper self energy 



E*(cj,k) via @ 



^^"''^^ = c.-k7(2m)-S-(..,k) ' ^'^^ 

where the spin indices have been suppressed. To 0{Cq) there are only two diagrams for 
the proper self energy, which are shown in Fig. |I](a) and (b). The filled circle represents 
the Cq interaction from Eq. (|l]). The Feynman rules for evaluating these diagrams can be 
found in Refs. [|], ^. Figure ^ shows schematically how the 0{Cq) contribution modifies 
the distribution. Since particles can be kicked out of the Fermi sea by the interaction, some 
occupation probability is moved to states above the Fermi surface. (Note that the second- 
order diagram shown is a Feynman diagram and so includes modifications of both particle 
and hole states.) 

If we use the equivalent Lagrangian Ca with a ^ 0, n{k) differs already at 0{a). The 
additional self-energy diagrams up to 0{aCo) are shown in Fig. |1] (c)-(h). Here the empty 
circle represents the induced three-body vertex oc aCo and the empty triangle the induced 
off-shell vertex proportional to a in Eq. (^). (The Feynman rules for these vertices can 
again be found in Refs. Q.) Equivalently, the occupation numbers can be calculated 
from energy diagrams with an operator insertion corresponding to Eq. (^. The appropriate 
insertion for a\^a^ on a fermion line with energy uj and momentum p is 

(-z)(27r)3 hm e^-^'^^^p _ ^^^^^ ^ (^4) 

where a and (3 are spin indices and the factor exp{iu!ri) ensures the correct ordering of 
operators. The corresponding diagrams up to 0{aCo) are shown in Fig. H(a), where we have 
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FIG. 3: Feynman diagrams for the occupation number n{k). (a) gives the contribution from aJ-Oj^ 
while (b) shows contributions from the additional term for a ^ 0. Note that diagram (a) (iii) 
represents only one of three possible insertions of a\^a^^. The square with a cross and the square 
with a circled cross denote insertions of Eqs. (U) and (|l6|), respectively. 



indicated the insertion of Eq. (0) by the square with the cross. (Note that the diagram in 
Fig. ^(a)(Mi) represents only one of three possible insertions of the operator aj.aj^.) Only 
the first two diagrams in Fig. ^(a) are nonvanishing. We find 



An(A;)(j) 
An(A;)(jj) 



Atiq: 

4:Tra Co 



k) 



'2^g{g - 1) 



A3 (27r)%-o^ 



lim / d^p / / d^qe 



x53(p _ \,)Go{v)Go{l)Go{p - q)Goil + q) . 



(15) 



This discrepancy is not surprising, since the definition in Eq. (|I2|) corresponds to the operator 
given in Eq. (0). Even the operator for the total particle number corresponding to the 
Lagrangian is different from the naive expectation. For the occupation number, however, 
there is no symmetry that can be used to construct the operator and so an ambiguity 
occurs. This simple exercise can be continued to higher order, with increasingly sophisticated 
diagrams. The same pattern recurs, with additional diagrams depending on a induced at 
each order, generating a dependence in n{k). This implies that n{k) is not an observable. 

We might take Eq. (|^) together with £ as a definition and transform the operator nk 
at the same time as C. For finite a this implies that one has to calculate the additional 
diagrams shown in Fig. ^b) where the square with the circled cross denotes an insertion of 
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where Uj, pj, and aj label the energy, momentum, and spin of of the line j, respectively. The 
occupation numbers defined this way are independent of a by construction. The first two di- 
agrams in Fig. ^(b) exactly cancel the contributions from the first two diagrams in Fig. |^(a) 
(cf. Eq (|T5|)), while the third diagram vanishes. However, we had no basis for the original 
definition, since EFT is model independent and makes no a priori assumptions on the dy- 
namics. Thus we conclude that occupation numbers (or even momentum distributions) 
cannot be uniquely defined in general. 
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FIG. 4: Schematic diagrams for the interaction of an external source J(x) [wavy line] with the 
fermion system, "knocking out" a particle. 



How does this conclusion fit in with the standard analysis of (e, e'p) experiments, where 
the cross sections measured are, by definition, observables? There are further complications 
in a finite system, where the momentum distribution differs from the occupation numbers, 
but we can illustrate the analogous situation in our model problem by introducing an external 
source J{x) coupled to the fermion number. [Thus J{x) plays the role of the Coulomb field 
of the virtual photon in {e,e'p).] If we were constructing an EFT of an underlying theory 
(such as QCD), we would expect to need the most general coupling consistent with the 
symmetries, but for simplicity we will assume that J{x) has nonzero coupling only to ip'^'ijj 
in the original representation C. 

The noninteracting cross section for a = corresponds to the diagram in Fig. §(a). The 
same cross section is obtained for a 7^ 0, but only if we include the contributions from the 
induced vertex to the final state interaction in Fig. ^(b) [and to the initial state interaction, 
which is not shown] and the vertex contribution from the modified operator in Fig. §(c). In 
general, there are always contributions of all three types, dressed with additional interactions 
order-by-order, which are mixed up under field redefinitions. Isolating Fig. ^(a) is a model- 
dependent procedure since it depends on a. 

Note that the ambiguities have a natural size, as discussed for an analogous shifting of 
contributions between two-body off-shell and three-body vertices in Ref. [^. An interesting 
question is whether the stark difference in occupation numbers between nonrelativistic and 
relativistic Brueckner calculations [^ can be explained by the ambiguity. Similar ambigu- 
ities occur in other areas of physics as well. In deep inelastic scattering, the physical cross 
section can be written as a convolution of quark and gluon distributions with coefficient 
functions determined by perturbative QCD. It is well known, however, that the distribu- 
tions and the coefficient functions are individually scheme and scale dependent. The scheme 
and scale dependence of auxiliary quantities such as the pion distribution in a nucleon was 
recently clarified in Ref. |jT6|. Another question of current interest is whether the condensate 
fraction in a Bose-Einstein condensate, which is essentially the occupation number of the 



condensate, can be measured |]17 |. 

The conclusion that the extraction of a momentum distribution from (e, e'p) cross sections 
is ambiguous because of final state interactions and vertex corrections (e.g., meson exchange 
currents) is not a surprise. While it is well known that such ambiguities are present, the 
usual assumption is that there is a "correct" answer that can be extracted from experiment. 
In a similar spirit, different ways of implementing the impulse approximation for the response 
of many-fermion systems have been analyzed in Ref. [|18]. In contrast, the EFT perspective 
clearly implies that ambiguities in the extraction of momentum distributions in (e, e'p) 
cannot be resolved by experiment. It is not only that the momentum distribution is difficult 
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to extract but that it cannot be isolated in principle within a calculational framework based 
on low-energy degrees of freedom. Rather, such auxiliary quantities can only be defined in 
a specific convention, like a particular form of the Hamiltonian, regularization scheme and 
so on. It can still be useful to discuss such quantities within a given convention. All true 
observables, however, can just as well be described in a different framework that adheres to 
different conventions. 
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